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■ A fluctuation theorem is proved for the macroscopic currents of a system in a nonequilibrium 
\ steady state, by using Schnakenberg network theory. The theorem can be applied, in particular, in 
, reaction systems where the affinities or thermodynamic forces are defined globally in terms of the 

cycles of the graph associated with the stochastic process describing the time evolution. 
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I. INTRODUCTION 



Out-of-equilibrium systems are ubiquitous in Nature and they play an essential role in many physical, chemical 
and biological phenomena QSISSIEIEHIIjS EH El E3 • In an isolated system, a nonequilibrium state will 



spontaneously relax toward the state of thermodynamic equilibrium. However, a nonequilibrium steady state can be 
maintained in open systems by exchanging energy or matter with thermostats or chemiostats. Such reservoirs generate 
thermodynamic forces - also called affinities |2( _ because they introduce inhomogeneities in the temperature, the 
pressure, or the chemical potentials of the different species of molecules. In turn, the thermodynamic forces or affinities 



generate fluxes of energy or matter across the system, which contribute to positive entropy production as described 



at the macroscopic level by nonequilibrium thermodynamics. 

The study of nonequilibrium steady states has a long history. Already in the XlXth century, nonequilibrium steady 
states were studied in linear electric circuits composed of resistors, capacitors, and inductors. In this context, Kirchhoff 
\ developed a network theory to determine the steady state solution E3- With the beginning of the XXth century, 
the interest shifted to electric circuits with nonlinear elements such as vacuum tubes and, moreover, fluctuations in 
electric currents such as the Nyquist- Johnson and shot noises were discovered, in particular, under the influence of 
Einstein's work on Brownian motion and its electric analogue |l4|. Since 1948, nonlinear elements as transistors and 
diodes have become dominant in electric circuits. These elements have macroscopic characteristic curves which are 
£f~) , not linear contrary to resistors. At the mesoscopic level, all the electric circuits present fluctuations in their currents. 
In this regard, the currents are random variables and the time evolution of the system is a stochastic process. The 
random variables are described in terms of a probability distribution. This latter obeys a master equation ruling the 
stochastic process. The random transitions are described in terms of the transition rates which may depend on the 
state variables. 

Similar considerations apply to chemical and biochemical reactions. In the sixt ies, Hill and coworkers developed 
stochastic models for the transport of ions across a membrane in biophysics [Tol ITa [T^ . At the beginning of the 
seventies, Nicolis and coworkers developed stochastic models of autocatalytic chemical reactions to investigate the 
effects of fluctuations on the macroscopic nonequilibrium instabilities Autocatalytic chemical reactions 

have reaction rates which depend nonlinearly on the chemical concentrations. However, we recover a linear description 
H ■ at the mesoscopic level since the master equation ruling the time evolution of the probability distribution is linear. 

The analogy between the nonequilibrium steady states of stochastic processes and linear electric circuits was used 
since the sixties by Hill [TEl ITfi Il7| , Oster and coworkers 0, 0] , and Schnakenberg . Schnakenberg developed 
a network theory using the methods by Kirchhoff for linear electric circuits in order to determine the steady state 
q , solution 0. In this network theory, a graph is associated with the stochastic process. The vertices of this graph are 
the mesoscopic states of the stochastic system and the edges the possible transitions between the states. This graph 
. ! can be decomposed into cycles. Schnakenberg made the fundamental observation that the ratio of the products of 
the transition rates along a cycle and its time reversal is independent of the mesoscopic states and only depend on 
the macroscopic thermodynamic forces or affinities maintaining the system out of equilibrium |lfjj . This observation 
allowed Schnakenberg to express several thermodynamic properties and, in particular, the macroscopic entropy pro- 
duction in terms of the cycles of network theory Later, the theory of these cycles has been further developed by 
the Qians IMEI. 

Recently, remarkable properties of the molecular fluctuations have been discovered in nonequilibrium systems. The 
time evolution of the molecular fluctuations can be characterized in terms of their large deviations with respect to their 
average. Several large-deviation dynamical relationships have been obtained such as the formulas of the escape-rate 



and hydrodynamic-mode theorie s 123L l24l I2M l2fl |27| . as well as the nonequilibrium work theorem [28l |2^ | and 
the fluctuation theorems 12(1 l3?l l38l l39l liol [III ] . The fluctuation theorems are symmetry 

relations originating from microreversibility for the probability distributions of fluctuating quantities, which are of 
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interest for nonequilibrium thermodynamics and transport properties. These new large-deviation relationships are 
under current investigations to explore the range of their applications. 

The purpose of the present paper is to apply the fluctuation theorem to the nonequilibrium stochastic processes 
described by Schnakenberg network theory. These processes include, in particular, models of electronic transport in 
mesoscopic conductors |42| , as well as chemical and biochemical reactions maintained out of equilibrium by differences 
of chemical potentials between the reactants and the products 0, IE IS 0, IS IU EH H> El • Chemical systems arc 
out of equilibrium if the concentrations of reactants and products are not in their equilibrium ratio fixed in terms of 
the standard free enthalpies of formation of the reactions, as described in chemical thermodynamics |2l. l8t 1 1 (| . The 
differences of free enthalpies between the reactants and products play the role of thermodynamic forces or affinities. 
In this regard, the thermodynamic forces have energetic and entropic contributions, and do not reduce to mechanical 
forces as for external electric forces acting on charged particles. Instead, the macroscopic thermodynamic forces or 
affinities may differ from the mesoscopic thermodynamic forces which depend on the mesoscopic states and are thus 
time-dependent fluctuating variables as the currents |sl lij. IT(ij]. In such systems, the macroscopic thermodynamic 
forces or affinities should be defined in terms of cycles, as shown by Schnakenberg [JJj . 

Here, we arc interested in the fluctuation theorem for the currents which are fluctuating variables at the mesoscopic 
level. Our aim is to give an algebraic proof of this theorem for the generating function of the currents associated with 
Schnakenberg's macroscopic affinities. For this aim, we have to identify Schnakenberg's cycles. This is carried out 
by considering the characteristic determinant of the operator giving the generating function as its leading eigenvalue. 
Our results can be applied to a large variety of processes including nonequilibrium chemical reactions and electronic 
transport in mesoscopic conductors |43l |44| . For non-reactive systems driven in their bulk by an external mechanical 
force or for non-reactive systems driven at their boundaries by time-independent chemical potentials, the fluctuation 
theorem for the currents already established in Refs. [33l l3?i l38| is recovered. 

The paper is organized as follows. We introduce the stochastic description in terms of the master equation in Section 
[n] where Schnakenberg network theory is summarized. We also illustrate the theory with an example of a chemical 
reaction network. In Section IlIII the generating function of the currents is introduced by using Schnakenberg's graph 
analysis and the fluctuation theorem for the currents is proved in this general framework. In Section llVI we discuss 
the connection between the fluctuation theorem for the currents and the fluctuation theorem for the action functional 
by Lebowitz and Spohn |37| . We also show how the fluctuation theorem for the currents can be used to measure the 
entropy production in some mesoscopic systems even if the microscopic transition rates are not known. Conclusions 
are drawn in Section Ivl 



II. MASTER EQUATION AND SCHNAKENBERG NETWORK THEORY 



A. Master equation and its basic properties 



Many nonequilibrium phenomena are successfully described at the mesoscopic level in terms of Markovian random 
processes. In some simple systems, such processes can even be exactly derived from the underlying deterministic 
dynamics |25l |26j by introducing an appropriate partition of the phase space. In other systems, such processes can 
be rigorously derived from the underlying Hamiltonian classical or quantum dynamics in some scaling limit |45| . 
Continuous-time random processes are ruled by an evolution equation, called the master equation, for the probability 
to find the system in a coarse-grained state u> at time t: 



= J2 [W>V)P(a/,i) - WM<o')P(<",t)] (1) 

p,LJ r 

The quantities W p (ui'\u>) denote the rates of the transitions lo'-^uj allowed by the elementary processes p = 1,2, r. 
Each one of these elementary processes may independently contribute to the entropy production so that it is important 
to separate them in the master equation. This is the case in nonequilibrium reaction systems where each reaction p 
contributes to the entropy production [S 0, |2(], ES OH • 

The master equation is known to obey a ff-theorem for the entropy 

S(t) = ]T S°(w)P(u>, t)-J2 P{", t) In t) (2) 

UJ UJ 

associated with the probability distribution P{u>,t) describing the state of the system at the time t [SEJ- S ^) 
denotes the entropy due to the statistical distribution of all the degrees of freedom which are not specified by the 
coarse-grained state u> |40j . For instance, if the coarse-grained state u> only specifies the numbers of the particles of 
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the different species, S°(oj) is the entropy of the statistical distribution of the positions and momenta of the particles. 
The second term is the contribution to entropy due to the statistical distribution P(u),t) of the numbers of particles 
[4fT |. The entropy is here calculated in the units of Boltzmann's constant fee — 1-38 10~ 23 J/K. The time derivative 
of the entropy is given by 

dS d c S diS . . 

It ~ ~df + ~dt ^ ' 

in terms of the entropy flux d e S/dt and the non- negative entropy production 

^ = \ E Jp(<»,<o')A p {u,<S)>0 (4) 

where 

J p {uj,lu') =P(u,t)W p (u\(j')-P(u}\t)Wp(u'\u) (5) 
is the current of the transition wAu' and 

, P(u),t)W p (w\w') 

a ^^ pWwM (6) 

the corresponding affinity 0, 0, |4j} ■ The 77-theorem asserts that the entropy production is always non- negative 
in agreement with the second law of thermodynamics and thus characterizes the irreversibility of the process. 
In a stationary state where dP/dt = 0, the master equation (JJJ can also be written as 



]Tj p (c» = (7) 

in terms of the currents JHJl, which is nothing else than the Kirchhoff current law |Tcj . 
In the equilibrium stationary state, the conditions of detailed balance 

P^{u)W p {u\u') = P eq (Lj')Wp{uj'\u;) (8) 

are satisfied for all the possible forward and backward transitions ui^cu . Consequently, the currents JSJl and the 
affinities (© as well as the entropy production (@J all vanish at equilibrium. 



B. Schnakenberg network theory 

Schnakenberg has shown that many of the fundamental properties of a nonequilibrium random process can be 
investigated and understood by carrying out the analysis of the graph associated with the master equation [To| . 
Moreover, Schnakenberg found that the nonequilibrium constraints of a system are related to the affinities of the 
cycles of the graph, obtaining in this way an early form of the fluctuation theorem |lf)j | . 

Let us first introduce the tools needed for the graph analysis of the random process. For a system ruled by the 

master equation (JIJ, a graph G is associated as follows £ 10j: each state u> of the system corresponds to a vertex 

p 

while the edges represent the different transitions ui^u allowed between the states. Accordingly, two states can be 
connected by several edges if several elementary processes p allow transitions between them. 

A method has been provided by Schnakenberg to identify all the cycles of a graph ^(J • This method is based on 
the definition of a maximal tree T(G) of the graph G, which should satisfy the following properties: 

(1) T{G) is a covering subgraph of G, i.e., T(G) contains all the vertices of G and all the edges of T(G) are edges 
of G; 

(2) T(G) is connected; 

(3) T(G) contains no circuit, i.e., no cyclic sequence of edges. 
In general a given graph G has several maximal trees T(G). 

The edges I of G which do not belong to T(G) are called the chords of T(G). If we add to T(G) one of its chords I, 
the resulting subgraph T(G) + I contains exactly one circuit, C/, which is obtained from T(G) + 1 by removing all the 
edges which are not part of the circuit. The set of circuits {Ci,C2, ...,C;, ...} is called a fundamental set. An arbitrary 
orientation can be assigned to each circuit C; to define a cycle G;. A maximal tree T(G) together with its associated 
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fundamental set of cycles {C%, G%, Ci, ...} provides a decomposition of the graph G. We notice that the maximal 
tree T(G) can be chosen arbitrarily because each cycle C; can be redefined by linear combinations of the cycles of the 
fundamental set. 

An orientation is given to each edge of the graph G. The directed edges are thus defined by 

e = u A J (9) 

Let / be a directed subgraph of G. The orientation of the subgraph / with respect to its edges {e} is described by 
introducing the quantity 

!+l if e and / are parallel, 
— 1 if e and / are antiparallel, (10) 
if e is not in /. 

where e and / are said to be parallel (resp. antiparallel) if / contains the edge e in its reference (resp. opposite) 
orientation. It will be useful for the following to write the affinities and the fluxes in terms of these numbers. In the 
noncquilibrium steady state, the affinity of the edge © is defined by 

P s t{u')W p {w>\uS) y ' 

The affinity of an arbitrary cycle C is then written as 

A{C) = Y J S e {C)A e (12) 

e 

where the sum is taken over all the edges e of the graph G. It can be equally written as 

A(C) = ^5 e (C)B e (13) 

e 

with 

Wp(uj'\uj) 

Indeed, the terms In [P st (oj) / P s t(u>')} exactly cancel each other when calculated on a closed path. Moreover, we can 
introduce a scalar product 

(C,Ct) = Si(C)Si{Ci) (15) 

where Ci is the graph of the fundamental set associated with the chord I as introduced above. It takes the value +1 
(— 1) if C contains the chord I with the same (opposite) orientation, and otherwise. The affinity of any cycle C can 
then be decomposed as 

A{C)=Y d {C,C l )A{G l ) (16) 
i 

where the sum extends over all the chords. This shows in particular that the affinity of an arbitrary cycle is a linear 
combination of the affinities of a fundamental set 10]. Finally, applying successively Kirchhoff current law (JJJ, the 
mean current traversing the edge e is written as: 



J e = Y J S e (Ci)J{Ci) (17) 



Nonequilibrium constraints are imposed to a system if the temperature, the pressure, or the chemical potentials 
differ between the reservoirs surrounding the system. These constraints are characterized by the global thermodynamic 
forces or affinities defined by the differences of the temperatures, pressures, or chemical potentials of the reservoirs. 
These global affinities do not directly appear in the transition rates W p {lo\lo') of the master equation £[J. These rates 
depend on the temperature, the pressure, or the chemical potentials of the reservoir responsible for the transition 
lu -A ui' and also on the states us and a/, although the global affinities {A-y} are defined by the differences of the 
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temperatures, pressures, or chemical potentials of the reservoirs and are independent of the particular states u> ov u>'. 
In this regard, the global affinities are macroscopic. 

As observed by Schnakenberg |lfj |. the global macroscopic affinities {A 7 } can be identified by considering the cycles 
{Ci} of the graph G associated with the stochastic process. The fact is that the ratio of the products of the transition 
rates along the two possible time directions of any cycle Ci of the graph is independent of the states composing the 
cycle and only depends on the global macroscopic affinities {-4 7 }: 



n 



W p {loW) 
W p (cj'\lj) 



(18) 



where e <E Ci denotes the edges © in the cycle C;. The conditions (|18|l are satisfied in a large class of processes 
including diffusion processes in lattice gases, nonequilibrium chemical reactions, and electronic transport in mesoscopic 
conductors \WL 1.371 143. |44|. We notice that there can exist more cycles C\ than currents 7 between the reservoirs. 
The reason is that the graph describes all the possible states and transitions at the mesoscopic level while the 
currents 7 between the reservoirs are typically macroscopic and fewer than the mesoscopic states. The affinities or 
thermodynamic forces A(Ci ) associated with the various cycles Ci of a graph G may thus take the same value for all 
cycles corresponding to the same current 7 between two given reservoirs: A(Ci) = A-y for all Ci G 7. The present 
paper is devoted to the class of systems where the thermodynamic forces or affinities are defined by the Schnakenberg 
conditions i|18|l . These conditions are weaker than in systems with external mechanical forces where the affinities can 
be directly identified at the level of the transition rates W p (lo\lo') themselves. Here, the affinities only appear after 
closing the cycles Ci, which significantly complicates the reasoning. 



C. Example: Schlogl's reaction scheme 



Our purpose in this section is to illustrate Schnakenberg's network theory with the example of the nonlinear chemical 
network 

(19) 
(20) 
(21) 

Particles enter into the system from three different reservoirs, A, B, and C. The particles of the species X are produced 
by three different reactions from so many reactants. The reaction constants are denoted by k p with p = 1,2, 3. Out 
of equilibrium, fluxes of matter will cross the system between the three reservoirs. In a stationary state only two of 
them will be independent and it is necessary to distin guis h them as they are independent macroscopic observables. 
According to the mass-action law |il IM 1(1 FtI VsL Wi TToL fTTl l40l l43| . the transition rates of these reactions are proportional 
to the concentrations and given by 
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where X is the number of particles of the species X, [•] 
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(22) 
(23) 
(24) 

(25) 
(26) 
(27) 



denotes the mean concentration in the reservoirs, and Q is the 
volume of the system. The transition rates are nonlinear functions of the internal state X. The state of the system 
is described in terms of the probability P(X, t) that the system contains the number X of particles of the species X 
at the time t. The master equation ruling the time evolution takes the form: 



dt 



P(X,t) 



EE 

p=l v=±\ 



\w p (X - v\X)P(X -v,t)- W P (X\X - v)P{X,t) 



(28) 
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We notice that neither the macroscopic currents nor the macroscopic affinities are apparent in the master equation 
contrary to the systems where the affinities are given in terms of mechanical forces. In order to identify them, we use 
the graph analysis of Schnakenberg 0. The graph of this stochastic process is depicted in Fig. ^ 




FIG. 1: Graph G associated with the reaction network ^19^210 . 
At equilibrium, the detailed balance conditions (jHJ hold which imply that 

(29) 



h[A] _ k 2 [B] _ k 3 [C] _ (X) cq 



k-l fc-2 fc-3 ^ 

and the equilibrium probability distribution P cq (A) is Poissonian of average {X) cq . Out of equilibrium, the detailed 
balance conditions are no longer satisfied and the probability distribution PtgiX) of the nonequilibrium steady state 
is known to be nonPoissonian for such nonlinear chemical networks 0,0, 0,0 111- F° r the present model, the steady 
state is given by: 

A T,nWn(x-l\x) 

p * m = p " (0) n s>;U-i! 

= p m s n fci[A]» + fc 2 [B](x - 1) + k 3 [C](x - l)(x - 2)0^ 

^ k-ix + k_ 2 x(x - l)^ 1 + k„ 3 x(x - l){x - 2)Q~ 2 1 ' 

An affinity can be associated with each reaction (p = 1, 2, 3) according to Eq. JBJ): 

= P st (X)W p (X\X + I) 

MX) - ln p st (x + m P (x + i\x) (31) 

These mesoscopic affinities depend on the internal state X, and thus fluctuates in time. Therefore, they do not corre- 
spond to the time-independent macroscopic affinities. These latter are instead given by the Schnakenberg conditions 
((TBj) . A possible maximal tree is given by choosing all edges corresponding to reaction p = 3. The chords correspond 
in this case to all edges associated with the reactions p = 1,2. The corresponding cycles are depicted in Fig. ^ They 
start from the state X and go to the state X + 1 by the edge p = 3 and return to the state X by the edges p = 1 or 
p — 2. The corresponding macroscopic affinities are given by 



W 3 (X\X + 1)W P {X + 1\X) 
W 3 (X + l\X)W p {X\X 



or 



- 41 = ln T — TTaT' ^2 = In r , (33) 

which are independent of the internal state X and thus constant in time as it should. 

This example shows that the mesoscopic affinities © depend in general on the state of the system and will thus 
fluctuate in time together with the corresponding currents (JSJ along a stochastic trajectory of the process. This 
problem is overcome by introducing the macroscopic affinities defined by the Schnakenberg conditions l|18|) . These 
latter no longer depend on the mesoscopic state and are thus independent of time as expected for nonequilibrium 
constraints from macroscopic reservoirs. 



III. THE FLUCTUATION THEOREM FOR THE CURRENTS 



We derive in this section a fluctuation theorem for the currents for the class of systems defined by the Schnakenberg 
conditions l|18|) . We give a different demonstration than in reference 01 . Here, we adopt a definition of the currents 
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based on a given maximal tree. This definition allows us more freedom than in reference [43J because the choice of the 
maximal tree can be used to derive several large-deviation relations in a given system. Moreover, we here consider 
general processes, not only chemical reactions. Finally, the present proof provides an efficient numerical procedure 
to compute the quantities of interest. The connection with the fluctuation theorem for the Lebowitz-Spohn action 
functional [37j is discussed at the end of this section. 

We now begin the derivation of the fluctuation theorem for the currents. The instantaneous current on the chord I 
is defined by 

+00 

3l (t) = Si(e n )6(t-t n ) (34) 

n— — 00 

where t n is the time of the random transition e„ during a path of the stochastic process. We use the convention 
that ji is oriented as the graph G since Si(e n ) is equal to (— )1 if the transition e„ is (anti)parallel to the chord I. In 
particular, the fundamental cycles are such that Si(Ci) = 1. This means that the cycles are oriented in the same way 
as the chords I but this is not restrictive. The current l|34[l is a fluctuating random variable. The so-called Helfand 
moment 0] associated with the current is defined as the integrated current by 

Gi(t)= f dt'ji{t') (35) 
Jo 

The quantity of interest is the generating function of the independent currents crossing the aforementioned chords 
of the system 

Q({Xi};{Ai}) = lim -l] n (fl-E^/.''«'»M) (36) 

t— >oo t 

where the sum is taken over all the chords and 

Ai = A(d) (37) 
We notice that the generating function (|36|l can be written as [24| 

Q({Xl};{A}) = lim -i m ( e -£< A ' G <«) (38) 

t— >oo t 

in terms of the Helfand moments (|35|1 . 

We shall first prove the fluctuation theorem for the currents (|34|l , according to which the generating function of the 
fluctuating quantities IpJST) obeys the symmetry 



Q({A,}; {A}) = Q({Ai - A;}; {At}) (39) 



under the Schnakenberg conditions i|18[l . 
In order to prove Eq. Q39[l . we define 



F t («,{A,}) = (e-Si^,(t) )w (40) 



as the mean value of exp [— A;G;(t)] with the initial condition that the system is in state lu at time t — 0. We 
want to derive an evolution equation for this quantity. The quantity F t +dt(y>, {A;}) can be calculated by considering 
all the possible transitions occurring in an infinitesimal time dt starting from state lu at time t — 0, before evolving 
the ensemble average for a time t starting from the state 10' reached at time dt. All these possible transitions must 
be weighted with their respective occurrence probabilities multiplied by the contributions of the exponential factors 
exp [— J2i A/G;(<)]. The quantity F t +dt{u, {A/}) can then be written as 



F t+dt (uj, {A/}) = ^Pdtiuj A uj')exp 



-y~]si(w -A lo')\i 



F t {J,{h}) (41) 



in terms of the orientation numbers (jlt)(l with a sum taken over all the chords. In Eq. (|41J) . the exponential factor 
thus comes from the contributions of the quantity exp [— XiGi(t)] during the transition lu —> to 1 . To simplify the 
notations, we introduce the quantities: 

exp (—A;) if uj — » lo' corresponds to the chord I in its reference orientation 
e -zpHw ) s g-Ei Si(io^tu )\i _ ^ expj^^) if ^ JL, u ' corresponds to the chord I in the opposite orientation (42) 

1 otherwise 
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which depend on whether a transition is taken along one of the chords {1} or not. Pdt{u — > w') corresponds to the 
probability to jump to the state u' in a time c?i along the transition p and starting from u). This conditional probability 
is given by 



P dt {uJ Au') = < 

Inserting this formula into Eq. 1)4 we get 

F t +dt{u, {A/}) = 



W P Hu/)A + 0(A 2 ) 



if a; 7^ a/ 



1-Epv ^ P (a;|a;')rft + 0(A 2 ) if 



W = W 



(43) 



£ W„(a;|w') e- 2 -^l"')F t (^,{A ; }) 



F t ( w ,{A / }) + 0(dt 2 ) 



(44) 



so that the time derivative of F t (u>, {A;}) is given by 
d 



^5;(w -A w')A; 



rf ^ t ( W ,{AJ) = ^ j W^lw'Jexp 

= L {M F t ( W ,{A,}) 
with the initial condition F t=0 (uj, {A;}) = 1. Whereupon we find that 

(e-E.Mo* *'*(*')) = ]Tp st ( w )F t ( w ,{AJ) = ^P st ( 



F t (w', {A,}) - W p (w|w')F t (w, {AJ) 



(45) 



(46) 



With our conditions on the transition rates, the Perron-Frobenius theorem is of application and there is a unique 
maximal eigenvector Vr^x (which is positive) 



L { x l} V {Xl} = -Q({Xi})V {Xl} 



(47) 



so that the limit in Eq. (|36|l exists and the leading eigenvalue of the operator (|45|) gives the generating function of 
the currents. 

We continue our analysis by looking directly at the eigenvalues given by the roots of the characteristic equation 
of the operator 1)45)1 . We will make the connection between the Schnakenberg graphical analysis and this equation, 
showing that the characteristic equation presents the symmetry 1)39)1 of the fluctuation theorem for the currents. 
Consequently, the eigenvalues of the operator ()45)) should also have the symmetry. However, the eigenvectors do 
not necessarily present the symmetry so that the long-time limit must be performed in order to obtain the leading 
eigenvalue and that the symmetry arises in the generating function according to the Perron-Frobenius theorem. 

The eigenvalues of an operator L acting on a finite state space are given by the roots of its characteristic determinant: 



det (l-sfj =0 
Here, we introduce the matrix M composed of the elements 



Mi. 



Li 

u 



s if i = j 

The characteristic determinant ()48)l is thus given by the following expression 



detM = Yl (-I)*™ M lp(1) M 2p(2) • • • M np{n) 

pfESym(n) 



(48) 



(49) 



(50) 



where the sum is taken over all the permutations of the symmetric group of n objects, Sym(n). This group contains 
n\ elements. 7r(p) is the parity of the permutation p and takes the value ±1 depending on whether the permutation 
can be decomposed into an even or odd number of transpositions. 
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Each term in the determinant l|50[) is thus a product of the coefficients of the linear equation 1|45[1 . A graphical 
representation of each term can be obtained by selecting for each vertex an edge (with the additional possibility to 
stay at this vertex). If an edge connecting two states is not present, the corresponding terms in the determinant will 
of course be absent. 

There is the additional constraint that there cannot be two edges pointing towards the same vertex. This results 
from the fact that the permutations form a group and are therefore inversible. This constraint implies that the 
determinant only contains terms corresponding to some cycles of the graph plus some fixed vertex (called full circuit 
in graph theory). A simple example is the case of a 3 x 3 matrix A where the determinant contains the term 011(122033 
corresponding to three fixed vertices, the terms 011023032, 013O22O31, 012021033 where there is one fixed vertex and a 
cycle of length 2, and finally the terms 012023031, 013021032 with no fixed vertex and a cycle of length 3. For general 
matrices of n x n, this property results from the representation of the permutations p £ Sym(n) in terms of the cyclic 
permutations of length 1, 2, n. 

We can thus write any term in the determinant as 

M aa ■ ■ ■ M hh (MijMjk ■ ■ ■ M ni ) ■ ■ ■ (M pq M qr ■ ■ ■ M zp ) (51) 

where there can be cyclic permutations of order 2, 3, ...,n. Now, for each term of this form, there exists its adjoint 
where all the cycles have been reversed 

M aa ■■■M hh {M l3 M 3k ■ ■ • M m ) R • • • (M pq M qr ■ ■ ■ M zp ) R = M aa ■ ■ ■ M hh (M in ■ ■ ■ M kj M 3i ) ■ ■ ■ {M pz ■ ■ ■ M rq M qp ) (52) 

Let us first consider the simplest case where there is only one type of transition between any two states, that is the 
variable p does not enter in the master equation. Any term My can then be written as 

M ij = L ij - S 6 ij = h^ W ^ ( k \- S iU = J (53) 

In this case, we can associate each term in the determinant with its adjoint and consider the quantity 

M aa ■ ■ ■ M hh (M.jMjk ■ ■ • M ni ) ■ ■ ■ (M pq M qr ■ ■ ■ M zp ) + M aa ■ ■ ■ M hh [MijMjk ■ ■ ■ M m ) R ■ ■ ■ [M pq M qr ■ ■ ■ M zp ) R (54) 

Introducing the notations W aa — — J^b W(u) a \ajb) and Wij = W(u>i\cjj), this term can be written as 

(W aa -s)---(W hh -s) (W ij W jk ---W ni )---(W pq W qr ---W zp )[l[( H e- z H-')) + Tj( e -4(c) J| e *(oV))] (55) 

C ljEC C ujGC 

where the product is taken over all the cycles {C} appearing in the decomposition Ij51|l . We also used that e -2 ^'"-* = 1, 
e -z(u\w ) _ e z(w \u) anc j e q ua ti on |j*H^ for the affinity of a cycle so that 

W ij W jk ---W ni A(c) (56) 
W in ---W kj Wji V ' 

Moreover the quantity Jl^ec e" 2 ^'" > can be written in terms of the scalar product (|15|) as 

■Q e -2(fcV) = e -T. l (c,c l )\ l ( 57 ) 

Now, we can introduce this expression along with the expression lltjfl for the affinity in Eq. I|55|l to obtain 

(W aa _ s ) . . . (w hh _ s ) (WijWjk ■ ■ ■ W ni ) ■ ■ ■ (W pq W qr ■ ■ ■ ^ p )[n^ E ' (C ' Cl)A ' + I]y E ' (c < CiMi e^ c < Ci > A < 

c c 

= [W aa - s) ■ ■ ■ (W hh - s) (Wi^k".^...^^^ < A '-*>] (58) 

hence the symmetry 

\t->Ai- Xi for all I (59) 



As it is valid for any term, the determinant also presents this symmetry. Consequently, the eigenvalues given by the 
roots of the characteristic equation will also be invariant under the transformation l|59|l . 
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It is important for the preceding reasoning to be valid that a term and its adjoint occur with the same sign in the 
expression of the determinant. This is guaranteed because a permutation and its adjoint have the same structure in 
cycles (i.e., the same number of cycles of length k, for all k) hence a basic property of the symmetric groups Sym(n) 
implies that they are in the same conjugation class so that they have the same parity. 

Now, we can consider the case where there is more than one mechanism of transition between two states. The 
reasoning is similar except that the terms Lij are now written as 



La — 



-Y<pY<k W p{u)i\u k ) if i = j 



11 XLpWpiu^e-^^) if i ^j 
A nontrivial cycle (LijLj k ■ ■ ■ L ni ) in the determinant is then written as 

(L ij L jk ---L ni ) = [j2W P M^)e- z ^ iM ] []T^ 2 ( W > fc )e-^>*> 



(60) 



[E^ ( w "M e ~ 



z Pr (oj n \uJi) 



pi 



EE ' ' 'E WpMl^e-^^WpMMe^^^ ■ ■ ■ WVKN)e- Zp " ( "" M (61) 



PI P2 



if the permutation is of order r. This is equivalent to summing over all the different possible cyclic paths V between 
the states (iOiiOj ■ ■ ■ u> n uji): 



(L %3 L, k ■ ■ ■ L ni ) = Y^WijWik ■ ■ ■ w m)v II e~ ZpHuJ ' } 

{V} uj,p<EC-p 



(62) 



where W%j stands for W Pl (uji\ujj) and the index V is a possible path between the states {u>iUjj ■ ■ -uj n uji) obtained by 
selecting one of the possible transitions p between these states: 



-r-, Pl P2 p3 
V =U!i > Ulj > ld-2 > 



• tO n ► LOi 



(63) 



Each of these possibilities also form a possible cycle hence we can make the same reasoning as before and associate 
a term with its adjoint. Its adjoint will also contain all possible paths between the states but in a reversed way so 
that the terms can be regrouped with each other to extract the affinity as in Eq. H58|) . Indeed, if there are v cyclic 
permutations in the decomposition H51|) . the quantity (|54|l can here be expanded as 



(W aa -s) ••■ (W hh -s) 

X E'"E { [(W ij W jk ---W ni ) Vl n e-^H-') 
{Pi} {P v } ^,pec Pl 

{W l3 W 3k ---W m )\ J] e 2 ^'^ 

u>,p£C-p 1 



{W pq W qr ---W zp ) Vv J] e~ ZpHul,) 
{W pq W qr ---W zp )\ [] e*>^')] } (64) 



uJ,p£C-p v 



with W aa = — J2 P J2b W p (w a \wb)- We can put in evidence the transition rates to get 

(W aa - s) ■ ■ ■ (W hh -s) E ' ' ' E {(^Wi* ■ ■ ■ W m ) Vl ■ ■ ■ (W pq W qr ■ ■ ■ W zp ) Vv 
{■Pi} {v v } 

x I J^J e - z p(^l w ') . . . J^J e -M<" l<"0 -)_ e --^(o-Pi) J^J e *p(w»K-) . . . e -^(Cj>„) J~J e z p (uii\uij) 



uj,pGC-p v 



ui,pGC-p 1 



uj,p€LC-p v 



} (65) 



where we used the same relations as before. Decomposing the quantities Y[ u pecv e ±Zp ^ i l" J ) and the affinities A(C-p) 
in term of the fundamental set, we get 

(W aa - s) ■ ■ ■ (W hh -s) E ' • ' E i W a W Jk ■ ■ ■ W ni)Pi ■ ■ ■ ( W Pi W ar ■ ■ ■ W zp ) Vv 
{Pi} {PA 

e -T.i{Cv 1 ,C l )\ l . . , e -Y. l (Cv k ,c l )x l +e y. l (c-p 1 ,c l )(\ l -A l ) . ..^dC-p^Ci^-A,) 

(66) 
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or equivalently 



(w aa -«)... (w hh -s) 

{•Pi} {Pt>} 



(WijW, fc • • • • • • (WyWjk ■ ■ ■ W ni ) Vv 



(67) 



The symmetry (I59f) is thus valid even in the more complicated case where there are several possible transitions 
between two states so that the corresponding characteristic equation also presents the symmetry (|59f) . The proof can 
be extended to Markovian processes on infinite state space either by using truncation or by replacing the characteristic 
determinant by (|48|l by a functional determinant, in which cases the symmetry is valid in the domain of convergence. 

Finally, we can regroup the different mesoscopic currents associated with the same macroscopic current j 7 and 
which correspond the same affinity Ay = A(C 7 ), as discussed previously. We thus set A; = A 7 for ( G 7 so that 



j 7 (*) = E and g t (*)=E g '(*) 

267 ZG7 



(68) 



We thus obtain our central result: 



Fluctuation theorem for the currents: Under the Schnakenberg conditions M8\) . the generating function 0/ the 
fluctuating integrated currents G 7 (t) obeys the symmetry 



Q({A 7 }; {Ay}) - lim -- ^ X - G ^) = Q({Ay - A 7 }; {Ay}) 

t— >oo t 

In the noncquilibrium steady state, the mean value of the current j 7 (t) is given by 



dQ 
<9A-, 



= lim 

{A 7 =0} £->oo t J 



dt'{j 1 (t')) = lim -<G 7 (t)> 



(69) 



(70) 



Its higher-order moments can be obtained by successive derivations, which shows that this function generates the 
statistical moments of the currents. We notice that reversing the affinities corresponds to reversing the sign of the 
currents and their fluctuations. 

The generating function l|69|) is related to a large-deviation property of the probability distribution of 
(1/t) J Q dt'j 1 {t') in the following way. The decay rate of the probability that the currents averaged over a time 
interval t take given values 



H({a 7 }) = lim -ilnProb 

t^oo t 

is the Legendre transform of the generating function: 



— J dt'jj(t') € (a 7 , a 7 + da 7 ) 



ff({a 7 }) = Max {AT> 



The symmetry relation l|69|l yields 



ff({-a 7 })-#({a 7 }) = ^A 



(71) 



(72) 



(73) 



or equivalently 



Prob 


{7 fodt'j-yit') S (a 7 ,a 7 + c?a 7 )| 




Prob 


{ 7 Jo dt'j 1 {t'') G (— a 7 , — a 7 + da 7 ) j 



exp AyQ! 7 i 



(t -> 00) 



(74) 



which is the usual form of a fluctuation theorem. 

We now have a fluctuation theorem for the macroscopic currents which are physically measurable quantities. We 
emphasize that, thanks to the above construction, this new symmetry in terms of the macroscopic affinities is com- 
pletely general in the range of validity of Schnakenberg conditions (|18|) . The present framework allows us to clearly 
identify the macroscopic currents between the reservoirs maintaining the system out of equilibrium. 
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Another important aspect of this construction concerns the numerical obtention of the generating function l|36|l . 
Indeed if we integrate Eq. ^45(1 . the quantity F t (oj,{Xi}) will decay in the long-time limit as 

X;^(w,{Ai})~e-«« A '»* (75) 

The generating function can thus be obtained from the rate of the exponential decrease of this quantity. Such an 
algorithm is important because it does not strongly depend on the nonequilibrium conditions. In contrast, a direct 
statistical approach requires a very large statistics in order to obtain the generating function. This statistics grows 
exponentially fast as the parameters A 7 differ from zero and as the nonequilibrium constraints increase. 



IV. CONNECTION WITH THE LEBOWITZ-SPOHN ACTION FUNCTIONAL 



In this section, we discuss the connection between the fluctuation theorem for the currents and the fluctuation 
theorem for the action functional by Lebowitz and Spohn |37j . These authors have introduced the action functional: 

z(t) = , W P1 (uj \uj 1 )Wp 2 (tJjVg) ■ ■ ■ W Pn (uj n -i\uj n ) 

measuring the breaking of detailed balance along a random path 

ui(t)=uiQ — ► u>\ — ► u>2 — ► > >jj n (77) 

The generating function of the fluctuating quantity (|76() is defined by 

q{ri) = lim — ln^^W) (78) 

i^oo t 

and presents the symmetry 

q{rj) = q(l - r,) (79) 
The mean value of the quantity 176fl grows at a rate given by the entropy production in a nonequilibrium steady state 

n Earns: 

^(0)=lim« = M (80 ) 

The connection between the two fluctuation theorems (|69f) and (|79f) can be established by the following reasoning. 
The action functional l|76|l can be written as 



Z(t) = Y J B e f dt'Ut') 
Jo 



in terms of the quantities 1)14(1. We here use the convention that the orientation of the fundamental cycles are such 
that Si(Ci) — 1. This means that the cycles are oriented in the same way as the chords I but this is not restrictive. 

If we associate a different parameter rj e with each possible transition, we can define the multi-variable generating 
function 

Q({77 e })= lim _iwe-S."« B «/o < ««'i«(*')) (82) 

J— >oo t 

We recover the generating function q{rj) by setting rj e — rj for all the edges e so that 

q(v) = Q(M) (83) 
The generating function l|82|) obeys the fluctuation theorem 

Q({Ve}) = Q({l-Ve}) (84) 
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as proved here below. We can absorb the local quantities B e to obtain the generating function of the currents on all 
the edges e as 

<2({A e }) = lim _i]n(e-£. *•/.?*''•(*')) = Q({X e /B e }) (85) 

t— >oo t 

which should obey the symmetry 

Q({X e }) = Q({B e - X e }) (86) 

Indeed, the evolution operator for this generating function can be derived as done in the previous section for the 
generating function of the independent currents. The only difference is that there is a different parameter A for each 
possible transition, and not only for the transitions corresponding to the chords I. This amounts to redefine the 
quantities z p (u)\u>') as 

e -* P M"') = exp [ - ^ S c {lo A u/)A e l (87) 

e 

where the sum is taken over all the edges e. Let us consider the case where there is only one type of transition between 
two states. The general case can be treated in a similar way. Using that 

a term in the characteristic determinant detM = det(L — si) plus its adjoint can be written as 

M aa ■ ■ ■ M hh (MijMjk ■ ■ ■ M ni ) ■ ■ ■ (M pq M qr ■ ■ ■ M zp ) + M aa ■ ■ ■ M hh {M tJ M jk ■ ■ ■ M m ) R ■ ■ ■ (M pq M qr ■ ■ ■ M zp ) K 

= (W aa - s) • • • (W hh - s) (W ij W j k---W ni )---(W pq W qr ---W zp )[]le-Ze s e^ (89) 

c c 

Using Eq. (|13fl . we can check that the transformation 

X e ^B e - A e (90) 

leaves this term invariant, which proves the fluctuation theorem (|86fl and thus l|84|l . 

We notice that, if the chords / are separated from the other edges e ^ {I}, the transformation 

A ( Ai - Xi 

X e -> -A e (e^{/» (91) 

also leaves invariant the expression l|89|l . The demonstration of the previous section can thus be applied and one finds 
that the generating function defined as 

Q({Xi}) = Q({X l },{X e = 0} e ^ {l} ) (92) 
obeys the fluctuation theorem l|39l) because 

Q({A*}) = Q({M,{A e = 0} e7 , w ) = Q{{Bi-X l },{B e } e ^ {l} ) 

= Q({Ai-Xi},{a}e*{i}) = Q({Ai-X l }) (93) 

where the second equality results from the symmetry (|9(J|I and the third from the symmetry (|91|l . Accordingly, the 
fluctuation theorem l)84|) or (|86|l implies 

Q({A 7 }) = QUA, - A 7 }) (94) 

which is the fluctuation theorem (|69|l for the macroscopic currents. This establishes the link between the generating 
functions and their corresponding fluctuation theorem. 

Finally, the entropy production can be obtained from the measurement of the probability distribution of the 
currents and using the corresponding fluctuation theorem. Indeed, one can first notice that the mean value of the 
action functional Z(t) can be expressed as the sum over all the mean currents crossing the chords multiplied by the 
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associated affinities. Precisely, thanks to equation i|13|l and our choice of orientation, we can write the affinities of the 
fundamental set as 



If we separate the chords from the other edges in the sum H81J1 and use Eq. Q95JI . we get 



Z(t) = J2 B i f dt'ji{t')+ Y, B - 1 dt 'Ut') 

I J ° e^{l} J ° 

= J2a f *'ii(0-EE^( c ')^ [ t dt' 3l (t')+ t B e [ dt'j e (t') 



e^{l} 



, Jo 



Using the property that 

Si(C v ) = Si, v 

if I and I' are chords, the term R(t) can be transformed as 



R(t) = -VV^(G)s e f dt'j l (t')+ V B e f dt'j e (t' 

r „/, Jo ,rn Jo 



e#{i} 



V B 8 [-V fdt'j l {t')+ fdt'Ut') 



(95) 



(96) 



(97) 



(98) 



Now, Kirchhoff current law l|17(l implies that the mean value of the rest (|98f) identically vanishes. Since the mean 
value of the action functional Z(t) increases with the rate of entropy production by Eq. I)8()[l. we hnd that the entropy 
production is given by the sum of the products of the macroscopic affinities with the currents: 



djS 
dt 



(99) 



By the fluctuation relation l|74|l . we can determine the affinities {A 7 } from the measurement of the probability 
distribution of the currents, which also gives the mean currents J7 7 . The entropy production can thus be obtained 
thanks to the formula l|99|l solely in terms of the nonequilibrium fluctuations of the currents. 



V. CONCLUSIONS 



In this paper, we have established the fluctuation theorem for the currents for stochastic processes in which the 
macroscopic thermodynamic forces or affinities are identified thanks to the conditions (|18fl of Schnakenberg ^(j ■ This 
class of stochastic processes includes models of electronic transport in mesoscopic conductors J4 4J4I , biophysical 
models of ion transport across membranes 0, 0, , and autocatalytic chemical reactions 0, IHlfiM Hi 01 HH • I n 
these stochastic processes, the mesoscopic affinities (©, which are often considered, depend on the mesoscopic state 
variables and are thus fluctuating in time, as illustrated for the case of SchlogPs reaction scheme in Subsection III CI 
Instead, Schnakenberg network theory provides a decomposition of the graph associated with the stochastic process 
in terms of cycles. The ratio of the products of the transition states along each cycle and its time reversal is directly 
related to a macroscopic thermodynamic force or affinity, as shown by Schnakenberg [Toj . This allows us to define 
the macroscopic affinities and the corresponding currents in the systems here considered. 

In the proof of the fluctuation theorem for the currents, Schnakenberg's cycles have been identified in the char- 
acteristic determinant of the operator giving the generating function of the currents as its leading eigenvalue. The 
cycles of the graph naturally appear in the structure of the determinant. In this way, we can take the Schnakenberg 
conditions 118|) into account and prove the fluctuation theorem for the currents. We have also proved that the fluc- 
tuation theorems for the currents and for the Lebowitz-Spohn action functional can be deduced from another one for 
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the generating function of the mesoscopic currents on each edge of the graph. This established a connection between 
the different fluctuation theorems, which concern different fluctuating quantities. Moreover, we have also shown that 
the fluctuation theorem and the measurement of the fluctuating currents should provide a new procedure to obtain 
the entropy production even if the transition rates are not known, as it is often the case in an experimental situation. 

In conclusion, this theory can be applied to a large variety of processes in physics, chemistry, or biology. The 
fluctuation theorem for the currents have important consequences for nonequilibrium nanosystems, which we hope to 
report on in a future publication. 
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